Incorporating molecular scale structure into the van der Waals theory of the 

liquid-vapor interface 
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We have developed a new and general theory of nonuniform fluids that naturally incorporates 
molecular scale information into the classical van der Waals theory of slowly varying interfaces. Here 
the theory is applied to the liquid-vapor interface of a Lennard- Jones fluid. The method combines 
a molecular field treatment of the effects of unbalanced attractive forces with a locally optimal use 
of linear response theory to approximate fluid structure by that of the associated (hard sphere like) 
reference fluid. Our approach avoids many of the conceptual problems that arise in the classical 
theory and shows why capillary wave effects are not included in the theory. The general theory and 
a simplified version gives results for the interface profile and surface tension for states with different 
temperatures and potential energy cutoffs that compare very favorably with simulation data. 



I. INTRODUCTION 

This paper applies our general theory of nonuniform 
fluids, described in several earlier publications jl], g, [| ^. 

to the liquid- vapor interface of the simple Lennard- 
Jones fluid. Our approach here can be viewed as a gen- 
eralization of the classical van der Waals (VDW) theory 
for the density profile of the liquid- vapor interface [?|] 
that (a) incorporates accurate thermodynamic data for 
the uniform fluid and (b) corrects the usual assumption 
that the interface profile is slowly varying. The new the- 
ory takes account of nonlocal molecular scale density cor- 
relations in a very natural way and can be applied to a 
wide variety of problems where the classical theory would 
fail. This perspective also provides a new and physically 
suggestive interpretation of the classical VDW theory for 
even a slowly- varying liquid- vapor interface that removes 
many of the conceptual problems and ambiguities that 
arise in standard descriptions 0. Thus it seems appro- 
priate to refer to it as a molecular scale van der Waals 
(MVDW) theory §. 

Since many aspects of the MVDW theory have been 
presented in some detail in previous work, here we will 
just outline the main features and focus on the new re- 
sults we find for the structure and thermodynamics of 
the liquid-vapor interface. In this specific application 
where the interface is often slowly varying, correction (a) 
plays the most important role and (b) is relatively less 
important, though conceptually (b) represents the most 
important advance and permits much more general appli- 
cation of the theory. In particular we will show how our 
theory, which first determines the structure of the liquid- 
vapor (LV) interface, can be used to calculate thermody- 
namic properties such as the surface tension. These re- 
sults will be compared to data from computer simulations 
i, §, In], O and to a simplified version of the MVDW 



theory that includes only correction (a) . We also discuss 
the role of capillary wave fluctuations (7| |l^, ^3), which 
should be taken into account when comparing theory and 
experiment. 



II. MOLECULAR FIELD APPROXIMATION 

An essential ingredient in the MVDW theory and in 
our interpretation of the classical VDW theory is the 
introduction of an effective single particle potential or 
"molecular field" to describe the locally averaged effects 
of the unbalanced attractive forces in the nonuniform 
LJ fluid |l4| . Since the attractive interactions are rela- 
tively slowly varying, such an averaged treatment seems 
physically reasonable. To that end, the LJ pair poten- 
tial w(r) = uq(t) + u±(r) is separated into rapidly and 
slowly varying parts associated with the intermolccular 
forces so that all the harshly repulsive forces arise from 
wo and all the attractive forces from u\. We describe the 
LV interface using a grand ensemble with a fixed chemi- 
cal potential ji l and temperature fc^T = giving two 
phase coexistence with bulk liquid and vapor densities p l 
and p v respectively in an external field <j)(r) = 0. 

The theory then approximates the structure of the 
nonuniform LJ system by that of a simpler nonuniform 
"reference" or "mimic" system at the same temperature 
but with a chemical potential /1q and purely repulsive pair 
interactions uo(r). These give the same repulsive inter- 
molecular forces as in the LJ fluid and are well approx- 
imated for most purposes by hard sphere interactions. 
The nonuniformity in the reference system is induced by 
an appropriately chosen effective reference field (ERF) 
4>r(t) that is supposed to take into account the locally 
averaged effects of the attractive interactions in the LJ 
fluid. 
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How should 4>R(r) be chosen? Since we want the refer- 
ence fluid structure to approximate that of the full fluid 
to the extent possible, it seems reasonable to determine 
4>r(t) formally by the requirement that the local (singlet) 
densities at every point r in the two fluids are equal fll5| : 

Po(r; [M> Mo) = P( r 5 (!) 

Of course this density is not known in advance, so in prac- 
tice we will make approximate choices for cf>R motivated 
by mean or molecular field ideas. Here the subscript 
denotes the reference fluid, the absence of a subscript 
the LJ fluid, and the notation [0] indicates that the cor- 
relation functions are functionals of the external field <fi 
(which in the present case is zero in the LJ system and 4>r 
in the reference system). Unless we want to emphasize 
this point, we will suppress this functional dependence, 
e.g., writing Eq. (Q) as p (r) = p(r). 

III. CLASSICAL VDW INTERFACE EQUATION 

As discussed in detail in [||, we can derive the classi- 
cal VDW interface equation from this starting point by 
making two additional approximations. We briefly dis- 
cuss this interpretation of the classical theory and then 
describe how our new MVDW theory improves on both 
approximations. 

A. Simple molecular field approximation 

First, the classical VDW theory uses the simple molec- 
ular field (MF) approximation for the ERF 4>r: 

^ii(ri) = 0(ri) + J dr 2 p o (r 2 ;[0fl],Mo) u i( r i2) + 2a//, 

(2) 

where 

a =-\ J dY 2 u l {r l2 ) (3) 

corresponds to the attractive interaction parameter a in 
the uniform fluid VDW equation, as discussed below. 
This is just a transcription of the usual molecular field 
equation for the Ising model to a continuum fluid with 
attractive interactions Mi (r) and can be derived in a num- 
ber of different ways jllj . The connection to the un- 
balanced attractive forces is perhaps most clearly seen in 
the derivation in 10,0], which starts from a formally ex- 
act description of the force balance in a nonuniform fluid 
and arrives at Eq. (0) by a series of physically motivated 
approximations. 

For the LV interface we have <fi(r) — 0, but it is conve- 
nient in what follows to keep a general <f> which we will 
then set to zero. In that case we will also choose Pq so 
that the density Pq of the uniform reference fluid with 



4>r = equals p . With this choice the ERFc^r vanishes 
on the liquid side far from the interface where the density 
becomes equal to p l . 

Another special case of Eq. (^) arises when <fi is a con- 
stant. Since a constant field in the grand ensemble is 
equivalent to a shift of the chemical potential, Eq. (j^) 
then relates the chemical potentials in the uniform LJ 
and reference fluids ||. Equation (0) thus yields the fa- 
miliar uniform fluid VDW result [|7J 

M(P) = Mo(p) - 2pa, (4) 

where p(p) and po(p) denote the chemical potential as 
a function of density p for the uniform LJ fluid and the 
reference fluid respectively. In the classical theory p(p) 
is defined for all p in terms of known reference system 
quantities by the right side of this equation. Since the 
uniform reference fluid is well defined for all densities be- 
low freezing, no problems arise from densities in the two 
phase region of the LJ fluid. The MVDW theory will use 
a slightly different expression for <pR in Eq. (|ll|) below 
that gives a more accurate description of the thermody- 
namics of the uniform LJ fluid. 

In general, to calculate 4>r a self-consistent solution 
of Eq. (|2|) is required, since 4>r appears explicitly on 
the left side and implicitly on the right side through 
Po(r; [4>r], Pq). Thus a useful implementation of the MF 
idea requires a way to accurately determine the density 
response po(r; [4>r]> Mo) induced by a given <pR. 



B. Local response to ERF 

The classical VDW interface equation results when 
a second approximation, appropriate only for a slowly 
varying field, is used to estimate the density response 
po(r; [</>_r], /Xq). This hydrostatic approximation for the 
density takes account only of the local value of the field 
through a shift in the chemical potential [|[ [|, ||] . Thus 
po(ri! [<Ar],Po) i s approximated for each ri by pJJ 1 , the 
local hydrostatic density response, which satisfies 

Mo(Po 1 ) = Mo - 0flOi)- (5) 

Hence the nonuniform density po(ri) at each ri is as- 
sumed to equal p^ 1 = po(i"i; [O^pg 1 ), the density of the 
uniform reference fluid in zero field at the shifted chemi- 
cal potential p^ 1 = pg — 0_r(i"i), given by the right side of 
Eq. (^). When <]>r{y\) is a constant, this gives the exact 
result. 

The superscript ri in p^ 1 is meant to remind us that 
Pq 1 , like p^ or pjj , represents the density of the uniform 
reference fluid at a particular chemical potential p^ 1 , 
which from Eq. (^|) depends parametrically on ri through 
the local value of the ERF. Thus when ri is in the bulk 
liquid (vapor) phase then pjj 1 reduces to Pq (pg). When 
4>r is very slowly varying this "local field" approximation 
is quite accurate and in this special case is equivalent to 
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the local density approximation made in the usual inter- 
pretation of the VDW theory J?], [n]. However this ap- 
proximation ignores the nonlocal excluded volume corre- 
lations that can be induced by a more rapidly varying 
This represents a major limitation of the classical theory 
in more general applications, and will be corrected in the 
MVDW theory. 



C. Classical interface equation 

The classical VDW interface equation follows immedi- 
ately when po( r 2) is replaced by in Eq. (|J) and the 
latter is substituted into Eq. [|. This yields an integral 
equation for Pq 1 , which from Eq. ([!]) is supposed to equal 
the density in the full LJ fluid: 



MoOS 1 ) = Mo - 0( r i) - J dr 2 p r 2 ui(r u ) - 2ap\ (6) 

This can be exactly rewritten in a more standard form 
using p(p) as denned in Eq. fl): 



MM) 1 ) =//-#•!) - / dr 2 [p 



Po>i(r 12 ). (7) 



Specializing to the case of the LV interface with planar 
symmetry and <j) — and expanding Pq 2 to second or- 
der in a Taylor series about pg 1 (consistent with the as- 
sumption of a slowly varying interface) yields the classical 
VDW differential equation for the interface profile pg: 



M(Po) -H = m 



dz 2 



where 



m — — / dvr u\(r). 
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(8) 



(9) 



Equations ([?]) and (|J) are completely equivalent to the 
VDW theory for the LV interface as it is usually pre- 
sented . In our derivation the theory describes hydro- 
static densities in the reference system, and p(p) is also 
defined in terms of reference system quantities given on 
the right side of Eq. (ji|). This provides a simple and con- 
sistent interpretation of the classical theory that avoids 
all the conceptual problems associated with densities in 
the two phase region of the LJ fluid that arise in tradi- 
tional approaches. 

In this derivation we have obtained the VDW interface 
equation directly, without first approximating the free 
energy. We will later show how to determine the interface 
free energy in this approach. First we will discuss our new 
MVDW theory for the interface profile, which improves 
on both approximations made in the classical theory. 



IV. MVDW THEORY FOR THE LV 
INTERFACE PROFILE 

In the limit of a uniform system, Eq. (||) describes all 
effects of attractive interactions in terms of the constant 



parameter a as in the van der Waals equation. While 
this very simple approximation captures much essential 
physics and gives a qualitative description of the uniform 
fluid thermodynamic properties it certainly is not quan- 
titatively accurate. In particular, when used to describe 
a slowly varying liquid- vapor interface, it will predict 
shifted (molecular field) values for the densities of the 
coexisting bulk liquid and vapor phases. The main prob- 
lem with the classical theory in this case is not so much 
its description of the local density gradients, which are 
often small, but its predictions for the thermodynamic 
properties of the coexisting bulk phases. The first cor- 
rection made in the MVDW theory is to modify Eq. (Q) 
so that it agrees with a given equation of state for the 
uniform system while still giving reasonable results for 
nonuniform systems S. 



A. Modified molecular field approximation 

To achieve quantitative agreement with known ther- 
modynamic properties of the uniform LJ system we can 
replace the constant a by a function a that depends 
(hopefully weakly, to the extent the van der Waals thcor 
is reasonably accurate) on temperature and density 
Thus, instead of using the MF approximation for p(p) as 
in Eq. (Q), we assume that p(p) is known from an accu- 
rate bulk equation of state. In particular, we determine 
p(p) from the 33-parameter equation of state for the LJ 
fluid given by Johnson, et al p8| . This provides a very 
good global description of the stable liquid and vapor 
phases in the LJ fluid and provides a smooth interpola- 
tion in between by using analytic fitting functions. Thus 
it naturally produces a modified "van der Waals loop" in 
the two phase region and seems quite appropriate for our 
use here in improving the simplest MF description of the 
uniform fluid. 

Now we relate this accurate p{p) to the known pa(p) 
through a function a(p) defined so that 



M(M) = Mo(p) ~ 2pa(p). 



(10) 



Thus the exact chemical potentials in the uniform LJ and 
reference systems are related in the same way as is pre- 
dicted by the simple MF approximation of Eq. (||) except 
that the constant a is replaced by a (temperature and 
density dependent) function a(p) chosen so that Eq. ( JTo| ) 
holds. We showed in || that the ratio a(p)/a is indeed 
of order unity and rather weakly dependent on density 
and temperature. 

Because of the strictly local response in Eq. (||) , these 
results for a constant field can also be used to determine 
exact results in the hydrostatic limit of a very slowly 
varying field. We want to modify Eq. (||) so that in 
the hydrostatic limit it will reproduce these exact values, 
while still giving reasonable MF results for more rapidly 
varying fields. 

There is no unique way to do this, but the following 
simple prescription seems very natural, and gives our fi- 
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nal result, which we have called the modified molecular 
field (MMF) approximation for the ERF [| : 

4>r(ti) - <f>(ri) = ' J dv 2Po(r2; [0j?J,Mo)"i( r i2) 

+MpV- (ii) 

Thus the molecular field integral in Eq. (0) is multiplied 
by a factor a(pQ 1 )/a of order unity that depends on ri 
through the dependence of the hydrostatic density Pq 1 
on the local value of the field <^r(ti), and the constant 
2ap is replaced by the appropriate limiting value of the 
modified integral. The MVDW theory assumes that the 
ERF is given by Eq. (|llj) rather than Eq. (§). 

B. Nonlocal response to the ERF 

The second important correction made in the MVDW 
theory is to determine more accurately the full nonlocal 
response po(ri] [4>r], Po) to tne ERF, thus correcting the 
local hydrostatic density p^ 1 used in the classical theory. 
We introduced in || a simple and generally very accu- 
rate method for calculating the structure and thermody- 
namics of the reference fluid in the presence of a general 
external field, using linear response theory in a locally 
optimal way to calculate the nonlocal corrections to the 
hydrostatic density. For a very slowly varying field the 
theory gives the hydrostatic density and for a hard core 
field the theory naturally reduces to the Percus-Yevick 
approximation |l9t| . 

The result is an integral equation for the density 
po(ri), which we refer to as the hydrostatic linear re- 
sponse (HLR) equation: 

po(ri) =p r o 1 +P1 1 Jdr 2 c (r 12 ;p^)[ Po (r 2 )~p^}. (12) 

Here co(ri 2 ; p^ 1 ) is the direct correlation of the uniform 
reference fluid at the hydrostatic density p^ 1 . This can 
be accurately approximated using known results for the 
uniform hard sphere fluid, as discussed in (|, |j. The ri 
dependence of the linear response kernel cq through p^ 1 
is the most important new feature of the HLR equation. 
A discussion of the ideas leading to Eq. (|l2|) and of its 
advantages over standard methods, along with numerical 
details of its solution, is given in [|| Ej . 

C. Two step method and the MVDW theory 

The MVDW theory for the LV interface arises from 
the self-consistent solution of Eqs. ( [TT|) and (|l2|). A 
two-step iterative method proved sufficient in all cases 
tested. Given a starting guess for (fo(ri) one computes 
in the first step the local hydrostatic density response 
Pq 1 from Eq. (||). Then in a second step nonlocal cor- 
rections leading to poiri) are determined from (|l2]), and 



this is used in Eq. ( pi] ) to give a new estimate for 0fl(ri). 
This process is iterated to self-consistency, and accurate 
numerical results are readily obtained. 

D. Simplified MVDW interface equation 

A simplified version of the MVDW theory arises when 
one skips the second step and assumes that po( r i) = Pq 1 
as in the classical theory, while still using the accurate 
equation of state to determine p(p) from Eq. ([lO]). While 
the full MVDW theory is straightforward to implement, 
the nonlocal corrections for the LV interface are often 
small and the use of the hydrostatic approximation allows 
for a more direct comparison with the classical th eory . 
Using this approximation and Eqs. (fTT|) , (||), and (|1C|) , 
we obtain an integral equation analogous to the classical 
equation (]?]). Expanding for simplicity to second order 
and assuming planar symmetry yields a generalization of 
the classical interface equation (||): 

p(p^)-p e = m a (p^)d 2 p^/dzl (13) 

where 

m a (p) = ma(p)/a. (14) 

Following Rowlinson and Widom Q(RW) we can de- 
fine 

-W(p) EE f(p)-^ P +p e , 

= p[p(p)- f /}-[p(p)-p% (15) 

Here f(p) is the Helmholtz free energy density from our 
analytic equation of state that corresponds to p(p) given 
by Eq. (|l0|), where p(p) = df(p)/dp and the associated 
pressure p(p) — —f(p) + pp(p) from standard thermo- 
dynamics. We see that W(p) vanishes in the coexisting 
bulk liquid and vapor phases and note that the left side of 
Eq. Q is given by -dWlp^ 1 ) / 'dp^ . By interpreting the 
latter as a "force" , Eq. ( |l3|) is analogous to Newton's law, 
with p^ 1 the "displacement" , z\ the "time" and m a (p) a 
density (or "displacement" ) dependent "mass" . 

Note that Eq. (|l^) differs from the analogous equation 
that would arise in the classical theory from assuming a 
density dependent "mass" in the gradient correction to 
the free energy. As shown by RW in their Eq. (3.10), the 
latter would generate an additional square gradient term 
in the interface equation (|l3|). It is difficult to see how 
such a term could arise naturally in our approach. 

We report results here for the even simpler theory that 
arises when m a (p) in Eq. ( |l3|) is replaced by its classical 
value m given by Eq. (^|). The resulting simplified in- 
terface equation has the same form as the classical equa- 
tion (||). However it uses the accurate expression for p(p) 
given byEq. @, which assures a proper thermodynamic 
description of the coexisting bulk phases, while (inconsis- 
tently) retaining the classical expression for m. This pro- 
duces some changes in shape of the interface profile when 



compared to that predicted by the MVDW theory, but 
preserves many qualitative features such as the depen- 
dence of the interface width on the thermodynamic state 
and the range of the attractive interactions. One major 
virtue of this approximation is that a very simple expres- 
sion for the profile p^ in terms of the inverse function 
z(po) (with arbitrary origin) follows immediately from 
Eq. ( |l3| ) by quadrature, as shown by RW: 

/rn \ 111 fPo 

z(p ) z(p v ) = (-) J ^ dp l-W(p)}- 1 / 2 . (16) 

Results from this simplified approach and the full 
MVDW theory will be discussed in Sec. |y] below. 

E. Nonlocal correlations and capillary waves 

The use of the HLR equation in the MVDW theory 
allows one to take account of nonlocal correlations in- 
duced by the ERF. These arise mainly from the packing 
of the harshly repulsive molecular cores and become sig- 
nificant at high density when the ERF is rapidly varying. 
In principle small amplitude excluded volume oscillations 
would be expected at high density from linear response 
theory far from any localized perturbation M, pl[ . The 
amplitude of these oscillations for a general liquid- vapor 
interface depends on the thermodynamic state and on the 
strength and range of the attractive interactions. These 
features control the steepness of the ERF (f>n(r), which 
mainly determines how significant these nonlocal correc- 
tions to the classical theory are in a given case. 

For our study here of the LV interface in the LJ fluid 
the attractive interactions are relatively slowly varying 
and these corrections are numerically small in most cases. 
However, near the triple point the interface is very steep 
and noticeable oscillations on the liquid side are predicted 
by the theory and, with much smaller amplitude, are also 
seen in the computer simulations. The classical local hy- 
drostatic approximation precludes a description of any 
such oscillations and gives no indication of where it can 
break down. So these corrections are conceptually im- 
portant even for "smooth" LV interfaces. 

The difference in oscillation amplitude between the- 
ory and simulation arises because the MVDW theory de- 
scribes a "static" interface in the reference fluid induced 
by the ERF. Any theory that takes account of attrac- 
tive forces only through an ERF and uses reference sys- 
tem correlation functions to approximate structure in the 
LJ system cannot properly describe the physics leading 
to the long-wavelength capillary wave fluctuations that 
occur at a real LV interface |7], E2[ These induce 

characteristic long-ranged pair correlations in the inter- 
face region of a real fluid that are completely different 
from the corresponding pair correlations in the reference 
system, which remain short-ranged for any reasonable 
choice of ^ii(r). In a sufficiently large system, the capil- 
lary wave fluctuations can wash out any excluded volume 
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FIG. 1: Density profiles of the LJ liquid- vapor interface 
for states with different temperatures and potential energy 
cutoffs. From left to right: T = 0.7 and r c = 5.0; T = 0.7 
and r c = 2.5; T = 0.85 and r c = 5.0; T = 0.85 and r c = 2.5; 
T = 1.1 and r c = 5.0. Symbols are results of Mecke et al. |]] 
Lines are predictions of the MVWD theory. 



oscillations at the real LV interface, and indeed the entire 
interface profile p(r) itself [Ql 

However the small system sizes studied in computer 
simulations or encountered in most experiments on inter- 
faces in confined geometries cut off most effects of such 
capillary wave fluctuations. It is reasonable to interpret 
the MVDW theory as providing one way of defining an 
"intrinsic" profile unbroadened by capillary wave effects 
and to compare its predictions directly to the simulation 
or experimental data, taking account of the residual fi- 
nite size capillary wave effects separately if necessary [[l| . 
This will be discussed further below. 

Workers using other approaches such as density func- 
tional theory or integral equation methods sometimes ar- 
gue that their theories describe correlations in the full 
nonuniform LJ system, and thus may include some or 
perhaps all of the effects of capillary wave fluctuations 
p2| , |23| . However such theories usually introduce approx- 
imations that relate correlation functions in the nonuni- 
form system to interpolated or weighted correlation func- 
tions in the uniform LJ system. One must then deal with 
the ambiguities arising from unstable uniform densities 
in the two phase region. We believe most such arbitrary 
schemes implicitly introduce a mean field character to 
the theory through the use of uniform fluid correlation 
functions that do not contain any capillary wave effects. 
However the precise physical implications of such approx- 
imations are very difficult to assess. In our interpretation 
of MF theory, both the strengths and the limitations aris- 
ing from the use of reference system correlation functions 
are clear from the outset. 
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V. RESULTS FOR INTERFACE STRUCTURE 
A. MVDW theory 

Figure 1 shows the interface profiles for the LJ fluid 
for states at three different temperatures. The dots give 
results of recent molecular dynamics (MD) simulations of 
Mecke et al. || . They made a careful study of the impor- 
tant changes in the interface profile and surface tension 
that arise from setting the force from the LJ potential 
to zero beyond a certain cutoff distance r c . Since there 
are unbalanced attractive forces in the interface region, 
effects from different cutoffs are much more important 
than in uniform systems, where the attractive forces es- 
sentially cancel. The first two curves on the left give 
profiles very near the triple point at a reduced temper- 
ature T of 0.7 with r c = 5 and 2.5 respectively. There 
are notable changes in the coexistence densities p v and 
p l from the different cutoffs. The next two curves give 
results at T = 0.85 with the same two cutoffs and the 
last curve on the right gives T = 1.1 with r c = 5. 

The lines give results of the MVDW theory, using the 
appropriate cutoff and shifted LJ potential. The effects 
of the cutoff on the bulk thermodynamics can be taken 
into account in the general equation of state of John- 
son, et al. jl|| used in the MVDW theory, and the ERF 
properly describes the averaged effects of the unbalanced 
attractive forces in the interfacial region. Overall there 
is very good agreement between simulations and the the- 
ory, which captures all qualitative effects of changes in 
temperature and cutoff radius. 

However, the theory predicts very noticeable density 
oscillations on the liquid side at the lowest temperature 
T = 0.7. As discussed earlier, the amplitude of these os- 
cillations is a sensitive function of temperature and cutoff 
and already by T = 0.85 their amplitude is greatly re- 
duced. Such nonlocal excluded volume correlations are 
to be expected when the density is high and the ERF is 
sufficiently rapidly varying. In other contexts they play 
an important role in the physics of nonuniform fluids . 

Indeed we find that the theory provide an exception- 
ally accurate description of po{r ; [</>_r], Hq), the density 
of the reference system in the presence of the ERF 4>r- 
This is illustrated in Fig. 2, which compares the theory 
(solid lines) for T = 0.7 and T = 0.75 directly to results 
from grand canonical Monte Carlo simulations we car- 
ried out for the reference system density (symbols) 
in the presence of the self-consistently determined ERF 
(dotted lines). Note that while the field is smooth, it is 
sufficiently rapidly varying in this case to produce a den- 
sity response with nonlocal excluded volume correlations. 
These are very well described theoretically by the HLR 
equation used in the MVDW theory and are completely 
missed by the local hydrostatic approximation (dashed 
lines) used in the classical theory. 

The differences in oscillation amplitude seen in Fig. 1 
arise because the reference system cannot describe the 
capillary wave fluctuations seen at a real LV interface, as 
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FIG. 2: MC simulation of the reference system in the pres- 
ence of the self-consistent field <f>R, given by the dotted lines 
(use the right vertical axis for fields) . The full lines and sym- 
bols denote densities as in Fig. 1. The dashed lines gives the 
hydrostatic density defined by Eq. (|^). The states are: (left) 
T = 0.7, r c = 5.0; (right) T = 0.75, r c = 2.5. 



discussed above. Some residual effects are present even 
for the relatively small system sizes used in the com- 
puter simulations. We can take them into account in an 
approximate way by using the standard prescription for 
Gaussian capillary wave smearing of an "intrinsic" in- 
terface. Thus we convolute the "intrinsic" interface as 
given by the MVDW theory in Fig. 1 with a Gaussian 
distribution of local interface positions h: 



P{h) = (2tts 



2s-l/2 



exp(-/i 2 /2s 2 )> 



where the width s of the Gaussian is given by 

2 k B T L 
s = In — . 

27T7 L s 



(17) 



(18) 



Here L is the lateral box size in the simulation and the L s 
is the short distance (wavelength) cutoff, which is sup- 
posed to be proportional to the bulk correlation length or 
the interface width j?], |l2[ |l3) . In this case the amount of 
smoothing depends mainly on the choice of L s , and the 
very reasonable choice of three times the interface width 
gives the data plotted in Fig. 3. The minimal choice of 
just the interface width gives too much smoothing when 
compared to the MD data. 

We see that the finite size capillary wave fluctuations 
have very little effect on the smooth profiles at higher 
temperatures, but are quite effective in damping out 
some of the "intrinsic" oscillations predicted near the 
triple point. Our purpose here is not to advocate this par- 
ticular and somewhat arbitrary prescription for smooth- 
ing the results of the MVDW theory, but to point out 
that while excluded volume oscillations are accurately 
described by the MVDW theory, their influence on the 
LV interface profile is reduced by capillary wave effects 
not captured by the theory. In many other applications, 
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FIG. 3: MVDW theory with capillary wave smoothing of 
profiles in Fig. f, with L s equal to three times the interface 
width. 



VI. SURFACE TENSION 
A. Basic formalism 

In the MVDW theory, we first determine fluid struc- 
ture. To calculate the interface free energy or sur- 
face tension 7 for the LV system, we proceed formally 
and imagine a special path where the density in the 
LJ fluid changes linearly p5| from that of the uniform 
liquid with density p l to the final LV interface profile 
p(r; [4> = 0],//) = p(r) as controlled by a coupling pa- 
rameter A with < A < 1: 



px(r) = p e + A[p(r) - p l 



(19) 



Here p\(r) = p(r;[(j>\], p e ) where 4>\(r) is the (gener- 
ally nonzero) external field that formally produces the 
partially coupled profile p\(r) defined by the right side 
of Eq. (]Tg|) . Since p\(r) = 5flx/S<))\(r), on integration 
the change in the Grand canonical free energy associated 
with this density change is exactly given by: 



e.g., structure near a hard wall or in a slit or pore, the 
capillary wave fluctuations are suppressed even further 
and the MVDW theory gives a good description of the 
excluded volume correlations such rapidly varying fields 
induce flal. 



B. Simplified MVDW theory 



Figure 4 gives the interface profiles predicted by the 
simplified MVDW theory from Eq. @. The theory 
seems to give very good results, perhaps even better than 
those of the full MVDW theory in Fig. 1! However part 
of this agreement is a result of a fortuitous cancellation 
of errors. As shown in Fig. 2, the hydrostatic approxima- 
tion used in the simplified theory completely suppresses 
the excluded volume correlations that should really be 
present at T = 0.7 in the reference system profile. As a 
result it seems to give better agreement with the simula- 
tion data for the LJ profile, where capillary wave effects 
also not properly described by the simplified theory wash 
out most effects of the oscillations. Note that the simpli- 
fied theory gives least accurate results at T — 1.1 where 
the interface is more slowly varying, and one would have 
expected the theory to be most accurate. As we will 
see later, results for the surface tension from the sim- 
plified theory are also less satisfactory. Nevertheless the 
simplified theory captures well the qualitative effects of 
temperature and cutoff on the interface profile and it is 
exceptionally easy to implement. 



A=l 



A=0 



dr 



dX p x {r) 



dX 



(20) 



Here — (3Q\ = In where 5^ is the Grand partition 
function for the system with field <f)\ . For the LV interface 
we have </>A=o( r ) = ^x=i(r) = 0. Since tt\ = o = —p e V, 
with p the coexistence pressure, which equals that in the 
vapor phase, the free energy difference on the left side is 
the desired interfacial free energy 7, and is independent 
of any choice of Gibbs dividing surface [Q. Integrating 
by parts, we have our basic starting point: 



»,r i'dXMr)^ 





dX 



or for the linear path: 



dr[p(r)-p*] / dA0 A (r) 



(21) 



(22) 



B. MVDW theory for surface tension 

The MF approximation used in the MVDW theory 
allows us to evaluate these expressions using reference 
system quantities. Thus we assume from Eq. ([!]) that 
P\(r) = Po\(r) = Po(r;[<f>Rx],Ho), where (/>r\(t) is the 
field in the reference system producing the same profile, 
formally related to 4>\ (r) by the MMF equation (ITT 



0ju(ri) = 0A(ri) + 
+2a(pV 



dr 2 Po\(r2)ui{r 12 ) 



(23) 



Using Eq. @ and Eq. (10), we can exactly rewrite 
Eq. ( p3| ) in a convenient form for use in Eqs. ( pl| ) and 





FIG. 4: Simplified MVDW density profiles from Eq. 
Same states as in Fig. 1 



(ri) = // 



PiP* 1 



OA i 



dr 2 [poA(r 2 ) 



Pol 



(16). FIG. 5: Surface tension 7 of the LJ liquid-vapor interface. 

The two sets correspond to the full LJ interaction potential 
(upper) and the cut-and-shifted LJ interaction with r c = 2.5 
(lower). Full lines are MWVD predictions. Dashed lines use 
the hydrostatic density to compute I n0 niocai from Eq. 
Dotted lines are results of the simplified MVDW theory from 
Eq. dio|). Results of simulations are from Mecke et al. j] 
(squares), Holcomb et al. Q (empty circle), Haye et al. |k] 

"1(7*12) -(filled circles) and Potoff et al. jjj (filled triangles). 

(24) 



In the MVDW theory p(p) is determined from the accu- 
rate equation of state and is given by Eq. (|l~0|), and all 
densities are calculated in the reference system. 

To calculate 4>\( r i) from Eq. (24) we start with the 
final self-consistent profile po( r ) and use Eq. ( |l9|) to de- 
fine po\ (r) . We then iterate the HLR equation (|12|) in an 
"inverse" way to find the hydrostatic density p \ associ- 
ated with a given p \ (ri ) (Although we know that po\ (r) 
is linear in A from Eq. (K%, unless the density is slowly 
varying this does not imply that the same condition holds 
for Pq X .) With this in hand, we use Eq. ( pi[ ) to determine 
</>a(>*i) for several intermediate values of A. The surface 
tension is then calculated by carrying out the integration 
in Eq. (E2h numerically. 



Simplified hydrostatic approximations for the 
surface tension 



Assuming that /9qa(i"2) — Pq\ in Eq. (|24|) we have 



a(p, 



OA / 



dr 2 [Pol -PoAW( r 12)> 

(25) 



and we will use this expression in Eqs. (|2lj) or (|22| ) 
to determine the surface tension. For A = 1 we have 
<^A=i(ri) = 0, and Eq. ( p5| ) reduces to the generalized 
interface equation ( |l3| ) on expanding p r ^ x to second order 
about Pq\ . The simplified MVDW interface equation dis- 

follows 
or m a by 



cussed earlier whose solution is given in E< 
on further approximating a(p) by a in Eq 
to in Eq. (|13|). 

Now let carry out the A integration in Eq. (|l]) using 
Eq. ( p5| ) for (f>\ (ri). A clear discussion of most of the tech- 
nical issues in given by RW. The first two terms on the 
right in Eq. (Eq) represent the local contribution to the 
interface free energy and can be integrated analytically: 




If we ignore the difference between poa(i"i) and p^ x as 
in the classical theory and assume the latter varies lin- 
early in A, then we require only the final profile Pn 1 and 
can carry out most of the A integration in Eq. ( |2l|) ana- 
lytically. We have already seen that this is a rather accu- 
rate approximation for the structure of the LV interface 
and we use this approximation only under an integral in 
computing the free energy. As we will see, this greatly 
simplifies the calculation of the free energy and also al- 
lows us to make contact with classical results and the 
simplified MVDW theory. 







dr / dX 



dr 



dp, 



dX 



[P -mOoa)] 



p* d Po\ 



OA' 



- drW(p r ), 



(26) 



where W(p) is given in Eq. (15). 

The last term in Eq. (Eq) gives the nonlocal contribu- 



tion to the free energy. Using Eq. (19) for p^ x we require 
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only pQ and from Eq. ( p2^ we find 

Inoniocai = J dr l [po ~ p l ] K { r l) J dv 2 W ~ Po' 1 M l ( r ^) 

(27) 

where 

K(r) = 2 [ dXXa(pl x )/a. (28) 
Jo 

Thus we have 

'f Ilocal ~\~ Inoniocai , (^^) 

and both terms require only p r Q . In the simplified M VDW 
theory discussed earlier, we set a = a or K{v)=l in 
Eq. (£7|) and expand p r 2 about p^ 1 to second order. Us- 
ing Eq. ( Jj"3| ) with to q = m, we see that these additional 
approximations imply that I non iocai = hocah as given in 
Eq. (|2(j). RW show in this case we can use the even 
simpler expression 

7 = / dp[-2mW(p)} 1/2 , (30) 
which does not require explicit knowledge of Pq. 

D. Results 

Figure 5 gives the surface tension predicted by the 
MVDW theory (with no capillary wave smoothing of 
the profile) , that arising from use of the hydrostatic ap- 
proximation as in Eq. (|2^) and (^), and that given 
by the simplified MVDW theory in Eq. @. We see 



that the MVDW theory gives very good agreement with 
the simulation data. The simplified MVDW theory is 
much less accurate. Since all theories give essentially the 
same results for Ii oca i, the problem with the simplified 
theory must arise from its treatment of nonlocal effects 
through the approximation Inoniocai = hocah which we 
see becomes increasingly inaccurate at lower tempera- 
tures. Equation ( p7|) provides a more accurate but still 
simple alternative to use of the full MVDW theory. 



VII. CONCLUSIONS 

The MVDW theory provides a simple and physically 
motivated approach that can describe the structure and 
thermodynamics of a fluid in a general external field. 
It optimally combines two standard approximations, a 
molecular field treatment of attractive interactions, mod- 
ified to give accurate thermodynamic data for the uni- 
form fluid, along with a linear response treatment of cor- 
relations in the reference fluid. The present application 
of the theory to the LV interface permits a new interpre- 
tation of the classical VDW theory that removes some 
ambiguities in standard treatments and shows how key 
features of the classical theory can be improved in a nat- 
ural way. The accuracy of the MVDW theory in this 
application provided additional support for the physical 
ideas behind the theory and for its quantitative utility. 

It is a pleasure to dedicate this paper to John Tully 
on the happy occasion of his 60th birthday. This work is 
supported by the National Science Foundation through 
Grant CHE-0111104. We are grateful to Yng-Gwei Chen 
and Jim Henderson for helpful discussions, and thank 
J. Winkelmann for sending us his computer simulation 
data. 
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